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We present some preliminary results for Bk at P — 6.0, using the twisted mass QCD formalism for the 
computation of bare matrix elements of the AS — 2 operator. The main advantage of the method is that mixing 
with other d = 6 operators under renormalisation is avoided. Moreover the operator renormalisation is performed 
in the Schrodinger functional (SF) framework, using earlier results of our collaboration for the corresponding step 
scaling function. 



1. Bk and tmQCD 

Let us start by recalling in brief the tmQCD 
framework for the computation of Bk 1 . We 
start from a continuum fermion action of the form 



+ s{x)[p + ms]s{x)} , (1) 



where — (u, d)'^ is a doublet of degenerate light 
quarks, and acts on isospin space. The axial 



transformation ip ^ e 



ia75T3/2 



leaves the form of the action invariant, and in- 
duces a rotational transformation of the mass pa- 
rameters. In particular, by choosing the rota- 
tion angle such that tan(a) = ^i/mi the stan- 
dard QCD action is recovered. The quantity 
Ml = mf + fi^ is left invariant, and can be iden- 
tified with the physical light quark mass. 

The above chiral rotation, on the other hand, 
can be seen as a change of fermion variables which 
induces a mapping between composite operators. 



*Based on a poster presented by P. Dimopoulos at the 
LATTICE 2003 Conference (Tsukuba, Japan). 



If we denote with a prime the quantities in stan- 
dard QCD variables then one has, in particular. 
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cos{a)Oy 



i sin(Q!)0: 



(4) 



where all the operators have an s — d flavour 
structure. These tree-level properties are imme- 
diately extended to the renormalised quantum 
theory, provided renormalised quark masses and 
composite operators are used in the above rela- 
tions. Then Eq. Q implies that, for the spe- 
cial case a = 7r/2, the physical matrix element of 
qAS=2 entering Bk is given in the twisted theory 
by a matrix element of the VA-I-AV part of the 
operator, which is protected from mixing under 
renormalisation by CVS symmetry and therefore 
renormalises multiplicatively 0. 

To compute Bk we regularise the theory using 
a Schrodinger Functional (SF) action with Wilson 
fermions. The action is nonperturbatively 0(a) 
improved in the bulk of the SF cylinder, and one- 
loop 0(a) improved at the time boundaries. The 



Set 


aMps 




Bk 


I 


0.3901(15) 


0.3892(16) 


1.022(11) 


II 


0.3561(20) 


0.3546(20) 


0.999(24) 


III 


0.3298(19) 


0.3284(19) 


0.984(23) 


IV 


0.3175(23) 


0.3040(24) 


1.003(27) 






Extrap. 


0.959(64) 



Table 1 

Values of the bare parameter Bk and pseu- 
doscalar masses for the four sets of quark masses 
and extrapolation of Bk to the physical kaon 
mass. (Sets correspond to 245, 150, 150, 102 
gauge configurations respectively.) 

bare value of Bk is extracted from a fit to the 
central plateau of the ratio of correlation func- 
tions 



(5) 



where ^ — [Za^o ~ iZyVf)] is the combina- 
tion of renormalised operators which represents 
the physical axial current in twisted variables, 
and O — L^^ ^ Cd(x)75Cs(y) is a pseudoscalar 
SF boundary source^. 

2. Results 

Here we present results at /3 = 6.0, all of them 
on a 16"^ X 48 lattice, corresponding to an ap- 
proximate physical volume of (1.5 fm)^ x 4.5 fm. 
All the simulations have been carried out on the 
APEmille machines at DESY-Zeuthen. Four sets 
of quark mass values have been considered. Three 
of them correspond to degenerate s — d quark 
masses ^, and the fourth one introduces a small 
(w 15%) s — c? isospin breaking. 

The lowest pseudoscalar meson mass can be ex- 
tracted either from an s — d two-point correlation 
function with the appropriate quantum numbers 
or, alternatively, from an s — d' correlation func- 
tion, where d' is an untwisted down quark with 
exactly the same mass as s. In Tabled (2nd and 

^In Eq. O, in standard SF fashion, the prime denotes the 
xo = T boundary. 

^The degeneracy holds for renormalised quark masses and 
up to 0(0^) cutoff effects. 



CQ 1.15 

DC 



is 



0.5 0.6 0.7 0.8 0.9 



Xo/T 

Figure 1. Rbk i^o) and the fit to the plateau for 
Bk for the data set I of Table ^ 

3rd columns) we show the values A/ps and Mpg 
for the corresponding pseudoscalar masses. If the 
twisted d quark is degenerate in mass with the 
s up to 0(0^) cutoff effects, then also the two 
meson masses should differ only by O(a^) effects. 
This is clearly seen in the data, which show no no- 
ticeable difference between A/ps and Mpg in the 
first three cases, whereas in the fourth case isospin 
breaking manifests itself in the non-equality of 
the masses. 

In Table n (4th column) we present our results 
for the bare value of Bk for each parameter set, 
as well as the result obtained by extrapolation 
to the physical kaon mass, which we take to be 
uMk — 0.2336. In all cases the fitting interval 
is Xo/a e [16,32]. As an example in Fig. ^ we 
depict the ratio Rbk ^nd the fit to the plateau 
which gives the Bk for the data set I. 

The data also allow to examine the chiral be- 
haviour of the {K^\0^^^'^\K°) matrix element, 
which can be directly computed from the ratio 



Bokk{xq) - 



{O'O) 



(6) 



and then studied as a function of {qMk)'^- The 
result is shown in Fig. |21 The operator O^f^ 
has been renormalised in the first of the SF 
schemes discussed in the next section, at the 
hadronic scale ^ = l/(2Lniax)- The intercept of 
the linear extrapolation is C = —0.00038(27), i.e. 
compatible with zero within 1.5 standard devi- 
ations. We interpret this as a signal of correct 
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Figure 2. Extrapolation to the chiral limit of the 
— K'^ matrix element. 

chiral behaviour, in view of the fact that the pre- 
diction of a vanishing matrix element in the chiral 
limit is only true for this quantity in the contin- 
uum. Here we are working with an unimproved 
matrix element computed at fixed lattice spacing, 
and the behaviour in Fig. [21 is still afi^ected by 
mild 0{a) effects. Moreover, as in previous stud- 
ies with Wilson fcrmions, the data lie in a region 
where LO chiral perturbation theory predictions 
have to be taken with care. 

3. Renormalisation 

For the nonperturbative renormalisation of 
Oya_+av '^^ techniques, as described in 

Nine diff'erent SF renormalisation schemes are 
available, each one of them supplying a formally 
(up to statistical correlations) independent deter- 
mination of the renormalisation group invariant 
parameter . Having only (3 = 6.0 results and 
thus no continuum extrapolation yet, an estimate 
of the latter is obtained by multiplying the bare 
value of Bk times the renormalisation constant at 
(3 = 6.0 and renormalisation scale l/2Li„ax, after 
which the result is multiplied by the continuum 
quantity B^°V-SF(l/(2imax)), obtained using 
the corresponding nonperturbative step scaling 
function and NLO perturbation theory (see Pj for 
details). The renormalisation constant and the 
ratio are scheme-dependent quantities, whereas 
the final result is not. In Table |21 we supply the 
resulting estimate for B^^ and B^^{2 GeV) in 



Scheme 


rRGI 


Bf^{2 GeV) 


1 


0.88(7) 


0.64(5) 


2a 


0.94(8) 


0.68(6) 


2b 


0.92(8) 


0.67(6) 


3a 


0.89(7) 


0.64(5) 


3b 


0.87(7) 


0.63(5) 


4a 


0.94(7) 


0.68(5) 


4b 


0.93(8) 


0.67(6) 


5a 


0.94(8) 


0.68(6) 


5b 


0.93(8) 


0.67(6) 



Table 2 

Values for B^^^ and i?^^(2 GeV) obtained after 
renormalisation aX (3 — 6.0. 

each of the nine schemes. 

The combination of the nine estimates of B^^ 
will ultimately yield a good control over system- 
atic effects and a reduction of the uncertainty of 
the final result. It has to be stressed, though, 
that the results in Tabled while being compat- 
ible within errors, involve different cutoff effects, 
and hence a reliable answer can be obtained only 
after extrapolation to the continuum limit (sim- 
ulations at higher values of (3 will follow soon). 
Note, however, that the (3 = 6.0 result is already 
close to the continuum limit values found in the 
literature, which suggests that cutoff results in 
our framework might turn out to be remarkably 
small. 
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